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Abstract 

At forward angles, and bombarding energies E > 200 MeV, the (p, n) and 
{n,p) reactions are thought to be directly proportional to the Gamow- Teller 



^^ I transition strengths in the nuclei. Assuming that this relationship also holds 

^D . for charge exchange induced by high-energy heavy ions, it would be very useful 

On ' in studies with radioactive beams. Contrary to this expectation, we show that 



the determination of Gamow- Teller and Fermi matrix elements from heavy-ion 
charge-exchange at forward angles is very inaccurate. 
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A major experimental and theoretical effort is under way to study the properties and 
applications of drip-line nuclei (for a review, see ref. [|I|). One of the main issues of this 
rapidly growing field is to understand the basic features of many unstable nuclei, which can 
only be assessed by using beams of radioactive nuclei. Of great relevance is the study of 
spin-isospin properties in unstable nuclei. A well-defined proportionality between heavy ion 
charge-exchange and beta-decay transition strengths would have important consequences. 
Such a relationship would be particularly useful to allow beta transition strength functions to 
be mapped out in drip-line nuclei. An accurate link between heavy ion exchange reactions 
and the well understood beta-decay processes would also put our understanding of the 
nuclear structure and reaction mechanisms of unstable nuclei on a surer quantitative footing. 

It is well known that (p, n) and (n, p) charge-exchange reactions at intermediate energies 
{E > 100 MeV) have been a powerful spectroscopic tool for measuring Gamow- Teller (GT) 
transition matrix elements @J^. This approach can also be used in reactions with radioactive 
beams, using a hydrogen gas target. However, due to the low luminosity of the beams, this 
process yields low counting rates. A more promising procedure would be the use of heavy 
ion charge exchange (HICEX) instead of {p, n) collisions. In ref. |Q] the possibility to use 
HICEX as a probe of Gamow- Teller matrix elements was suggested and a relation between 
the forward cross section and the Gamow- Teller matrix elements was found. A possible 
drawback of this approach is the complication originated by the internal structure of the 
target nucleus. This method was used as a guidance in ref. to draw conclusions on 
Gamow- Teller (and Fermi) matrix elements in ^^C and ^^A^ from the mirror nuclei reaction 

13^(13jy^ 13^)13^_ 

In this article we develop a strong absorption model of HICEX based on a single particle 
approach. We use this simple structure model for qualitative purposes only. Our aim 
is to access the reliability of the relationship of HICEX measurements and the Gamow- 
Teller and Fermi matrix elements. The key point to be discussed here is that if we cannot 
show the existence of a valid relationship for such a simple model and for rather simple 
reaction partners, then very probably it would not exist for more complicated nuclear models. 
Because of the use of mirror nuclei, and their isospin symmetry, the experimental data on 
the ^^C{^^N, ^^Cy^N reaction is also an ideal test of the theory. 

In heavy-ion charge exchange reactions it is important to know the underlying nature 
of the exchange mechanism. In ref. 0] it was shown that in heavy ion collisions at high 
energies {E > 100 MeV.A) the exchange mechanism is dominated by the fundamental tt- 
and p- exchange. The differential cross section for the process is given by 
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where fi is the reduced mass of the system (1+2), k {k') is the initial (final) center-of-mass 
momentum, Ji and J2 are the angular momenta in the entrance channel, and the sum is an 
average over initial angular momentum projections and a sum over final angular momentum 
projections, denoted respectively by ai and 02- 
The transition matrix Tjj is given by 



T^f = 7^ld'R X^-^*m X^^\R)}ld'q e'""-^ Mia,, «2, q) 



(0.2) 
is the 



where R = (h, z) is the relative coordinate between the nuclei, x*' H-f^) X^'^H-R- 
incoming (outgoing) scattering wave obtained in terms of the optical potential for elastic 
scattering, and 



M = Jd'n dh2 <^Pi(ri) 5p2{j^) e^'i-(^^-^^) K.c/.(q) , 



(0.3) 



where 5 pi is the transition density of the nucleus i, which includes spin, isospin and spa- 
tial transitions. The meson-exchange potential Vqbep is a sum of vr-, p-, and w-exchange 
interaction. We use a modified meson-exchange interaction to include a non-spin-fiip term: 

Vexch = VoBEP + VrTi.T2 [0,|]. 

The eq. p.3| is separable in terms of the spin and isospin variables, yielding M. = 
M12 + M2^i, where 
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In the equation above the sums over /i and p! run over the spherical components of the 
spin vector. Fo-^(q), G'o-^(q) and -f^r(q) are functions of the central (L=0) and the tensorial 
(L=2) component of the charge-exchange potential. 

Since the initial and final wavefunctions are decomposed into their spatial and spin- 
dependent parts, it is appropriate to expand the exponential functions in eq. ^^ into multi- 
poles and use the Wigner-Eckart theorem, to obtain sums over products of Clebsh-Gordan 
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coefficients and reduced matrix elements of the form < <pji \\ii{qr)[a ® F/]// t^ '\\(t)j >, 
< 4>p\\Ji{Q^)['^®^i]r\\4>j >i and < (Jjji \\ji{qr) T^^^Wcp^ >. These matrix elements can be 



calculated using eqs. (A. 2. 19) and (A. 2. 24) of ref. §]. They involve geometric coefficients 
and overlap integrals of the form 

J,,y (g) = J R,,{r) R,{r) jj{qr) r^ dr , (0.5) 

in terms of the single-particle wavef unctions Rj{r). The cross section results from interfer- 
ence of terms corresponding to several possibilities of angular momentum transfer. 



From |0.4| it is straightforward to show the relationship between the heavy-ion charge- 
exchange cross sections and the Gamow- Teller and Fermi matrix elements. In beta-decay 
experiments a similar matrix element appears, but with the exponential factors replaced by 
exp 2(pe + Pi/).r, where Pe and p^ are the electron and the neutrino momentum, respectively. 
This exponential factor can be approximated by the unity, if |r| < Ti/pe and h/p,j. Since 
the energies of the light particles involved in beta-decay are of the order of a few MeV, the 
Compton wavelengths of the electron and the neutrino are approximately equal to 10'^ fm, 
i.e., about hundred times larger than the nuclear sizes over which the integration occurs. This 
contrasts with charge exchange experiments, where the exponential factor in ^^ involves the 
Compton wavelength of vr and p. Thus, the relation q.r ^ 1 is not well justified for HICEX. 

Assuming that the "low-g" (or long-wavelength) approximation is also valid for HICEX, 



i.e., inserting eq. p^ , with the exponentials replaced by the unity in eqs. |0?T| - p.4| , one gets 
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B{GT) = — Y. |<aVTW|a> and B{F) = — ^ | < a'|r(±)|a > (0-7) 

•J a' ,a J a', a 

are the Gamow- Teller and Fermi strengths, respectively. The coefficients Ci, C2, and c^ 
are products of geometric factors and integrals over the q- and /^-coordinates. The "low- 
q" approximation ^^ also includes a sum over all single-particle transitions, which is not 
explicitly shown for simplicity. 

Eq. p.6| is the result that one would like to have. Assuming that the charge-exchange 
potential (corrected for medium effects, etc.) is accurately known and that the same is true 
for the scattering waves, then the coefficients Cj could be reliably calculated. An experimental 
measurement of HICEX would then provide the magnitude of B{GT)^s and S(F)'s for 
several nuclei. This would be a specially important technique for studies of dripline nuclei. 
We now investigate if this relation is valid using the ^^C{^^N, ^^GY^N reaction at 105 
MeV/nucleon ||^ for comparison. 



Following ref. we normalize the interaction so that at zero momentum transfer 
VoBEP = Var = 265 MeV.fm^ and Vr = 100 MeV.fm^, appropriate for proton laboratory- 
energies of order of 100 MeV. We describe the scattering waves in terms of eikonal wavefunc- 
tions. Since an optical potential for this reaction at 100 MeV/nucleon is not known, we use an 
optical potential which fits the reaction ^^C+^^C at 85 MeV/nucleon [|ll|;0- To calculate the 
single-particle wavefunctions we use a potential of the form V = Vof (r) + Vso(J-0') {df /dr)/r, 
with /(r) = 1 + exp (r — Ro)/a , for the neutrons. For the protons we add the Coulomb 
potential of a uniformly charged sphere of radius Rq. For both ^^C and ^'^A^ we choose the 
set of parameters Vq = —53 MeV, Rq = 2.86 fm, a = 0.65 fm, and Vso = —15.5 MeV. In this 
simple model, the relevant transitions are lp^/2 — ^ ^Pi/2 (^ ~ ""-5 ^^ Py i-^-5 ^^^ proton in 
13 jY goes to a neutron in ^^C, and vice- versa) involving the operator r (Fermi transitions). 
For Gamow- Teller transitions (ar operator) the relevant transitions are Ip^^ — ^ ^Pi/2 ^^^ 

lPl/2 ^ 1^3/2- 

In table 1 we show the matrix elements of the single particle transitions: -B1/2- = 
Bi/2- — > Bi/2- and -B3/2- = -B1/2- — ^ -B3/2- for Fermi and Gamow- Teller transitions. 
For comparison, we also give the values of these matrix elements, in the case that the over- 
lap integral includes a spherical Bessel function of zeroth-order, joiqr), for additional two 
values of q. The matrix elements are drastically reduced for large values of q, especially for 
lPi/2 — ^ 1^3/2 transitions. 

In figure 1 we compare the experimental |^ and the theoretical calculation of the dif- 
ferential charge exchange cross section for the double excitation case, i.e. when both target 
(^^C) and projectile (^^A^) go to the 3/2~ state. Both experimental data and the calculation 
are given in arbitrary units. Thus, only the angular dependence of the cross section can 
be discussed. The angular dependence is rather insensitive to the structure model for the 
nuclei, but strongly sensitive on the parameters of the optical potential used to calculate the 
scattering wavefunctions. Only when such parameters are known from, e.g., elastic scatter- 
ing one can have a better knowledge of the influence of the Gamow- Teller and Fermi matrix 
elements on the charge exchange cross section. This, of course, assuming that the exchange 
potential is well known. In fact, the form of the exchange potential is crucial for this aim. 

In figure 2 we plot the ratio between the exact calculation and the low-q approximation, 



eq. |0.6| . With the single particle basis the expansion can be stopped at the / = 5 term. In 
fact, the term with 1 = dominates. This is shown in the figure where the contribution of 
the / = (dashed curve) term of the expansion is shown. The total sum includes terms up 
to / = 5 (solid curve). However, note that the ratio between the exact calculation and the 
low-q approximation amounts to a factor of 5 at the forward direction. This difference can be 



understood in a simple way by using eq. pl2| . If strong absorption and final state interactions 
could be neglected, x'^^^*(R) x'^+^(R) ~ exp jiQ.RJ-, where Q = k' — k is the momentum 
transfer to the cm. scattering. Thus, the integral over R yields a delta function, (5(q — Q). 
That is, Q = (forward scattering) would imply q = 0, too. Very forward scattering would 
then justify replacing the exponentials in eq. 2]4| by unity, and an equation of the form 
of eq. |0]^ would be obtained. However, strong absorption modifies this scenario: Q is not 
directly related to q, and the equation p.6| is a bad approximation. The larger the absorption 
radius, the worse this approximation should be. This may somewhat explain why for {p, n) 
reactions the low-q approximation is better justified, since the absorption radius is smaller 
in that case. But, even for (p, n) reactions the low-q approximation seems to be ill-defined, 
as was shown recently |]T^, especially for the weak transitions. 

Another possibility to get direct information on the Fermi and Gamow- Teller strengths 
from HICEX reactions would be in the case that the expansion of the exponential function 
in p.4| can be stopped at the term with / = 0. One gets matrix elements of the form 



< 0J'7'm'I^A'^''io(gr)|0J,„ >= j(jlm/i|jW)JF,y(g) 
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and a similar expression for Fermi transitions. The function J-'{q)jj> = \Iojj'{q)/Iojj'{q = 0)p 
has to be calculated with accuracy in this case. Similar arguments are given in ref. for 
{p, n) reactions. 

In figure 3 we show the percentage deviation between the exact calculation and a calcu- 
lation done with only the / = term of the expansion. Also shown (dashed curve) is the 
differential cross section in arbitrary units. At 6 = the difference is about 2% and increases 
largely close to the minima of the cross section. This is because, close to the minima (which 
are basically determined by the Iq term) the spherical Bessel functions j/>o contribute more 
to the integrals. 

Although the (J = 0)-approximation is not bad at 6 = compared to the full calculation, 
it cannot be considered as a useful tool to obtain the Fermi and Gamow- Teller strengths 
from HICEX, unless a good model for the overlap integrals can be done. This is however just 
what one would like to extract from the experiments, and not from theory. As one sees from 
table 1 the matrix elements involved in HICEX are strongly modified by the short range 
property of the exchange interaction. Only a small part of the transition densities, close to 
the surface of the nucleus is relevant for the matrix elements. In contrast, in beta-decay 
the whole transition density is relevant. To show this we decrease artificially the masses 
of the pion and the rho by a factor a. These masses enter into the form of the exchange 
interaction and are responsible for its range [^,0]. As we see from figure 4, as a increases, 

6 



increasing the range of the interaction, the ratio between the exact calculation and the low-q 



approximation decreases, reaching unity at a ~ 10. This means that eq. |0.6| would be valid 
for an interaction of range Tq ~ 10 x (/i/m^c) ~ 10 fm. 

In conclusion, we have shown that it is very difficult to obtain reliable estimates of 
Gamow- Teller or of Fermi strength from HICEX experiments. This might be surprising, 
since weak interaction strengths determined from {p, n) reactions are astonishingly good, 
although they have been shown to be inaccurate for important Gamow- Teller transitions 



whose strengths are a small fraction of the sum rule limit 1 12 
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TABLES 
TABLE L Modified Fermi and Gamow- Teller strengths for single particle transitions and as a 
function of q. The Gamow- Teller and Fermi operators are replaced by their respective product with 
a spherical Bessel function jo{qr). For q = 0, Jo{qr) = 1, and the usual Fermi and Gamow- Teller 
matrix elements are restored. 
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Figure Captions 



Fig. 1 - Differential cross section of tlie cliarge excliange reaction ^^C(^^A^/'^ C)^^C at 
105 MeV/nucleon for the double excitation case (target and projectile go to tlie 3/2~ state). 
Experimental data [^] and calculation are given in arbitrary units. 

Fig. 2 - Ratio between the exact calculation and the approximation p.6| for the charge 
exchange reaction ^^C(^'^A^/^ C)^^C at 105 MeV/nucleon, as a function of the scattering 
angle. 

Fig. 3 - Percentage deviation from the exact calculation and the approximation obtained by 
keeping only the / = term of the multipole expansion p.4| for the charge exchange reaction 
13(7(13 jY^is (jy^(j a^ X05 MeV/nucleon, as a function of the scattering angle. 



Fig. 4 - Ratio between the exact calculation and the approximation ^^ as a function of 
an scaling factor a for the charge exchange reaction ^^C(^^A^,^^ C)^^C at 105 MeV/nucleon. 
The scaling factor a is used to reduce the pion and the rho masses in the charge exchange 
interaction. 



